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We investigate effects of multipole moments in /-electron systems both from phenomeno- 
logical and microscopic viewpoints. First, we discuss significant effects of octupole moment on 
the magnetic susceptibility in a paramagnetic phase. It is found that even within mean-field 
approximation, the magnetic susceptibility deviates from the Curie- Weiss law due to inter- 
actions between dipole and octupole moments. Next, we proceed to a microscopic theory for 
multipole ordering on the basis of a j-j coupling scheme. After brief explanation of a method 
to derive multipole interactions from the /-electron model, we discuss several multipole or- 
dered phases depending on lattice structure. Finally, we show our new development of the 
microscopic approach to the evaluation of multipole response functions. We apply fluctuation 
exchange approximation to the /-electron model, and evaluate multipole response functions. 

KEYWORDS: octupole moment, j-j coupling scheme, Curie-Weiss law, fluctuation exchange approxima- 
tion 



1. Introduction 

Novel magnetism and exotic superconductivity in rare- 
earth 1 and actinide 2 compounds have attracted much at- 
tention in the research field of condensed matter physics. 
In order to elucidate the mechanism of these phenom- 
ena, it is quite important to describe /-electron states in 
a systematic way by considering correctly the symmetry 
of /-electron orbitals. On the other hand, we encounter 
difficulties to express /-electron states by means of a sim- 
ple combination of spin and orbital degrees of freedom, as 
has been done for d-electron states, since the spin-orbit 
interaction is strong for / orbitals in general. For such a 
case, it is rather convenient to exploit multipole degrees 
of freedom, including charge (rank-0), dipole (rank-1), 
quadrupole (rank-2), and octupole (rank-3) components. 
Fortunately, among such multipole moments, dipole and 
quadrupole correspond to ordinary spin and orbital de- 
grees of freedom, respectively. Thus, phenomena con- 
cerning dipole and/or quadrupole moments are basically 
understood by analogy with (^-electron physics. In this 
sense, octupole is considered to be the first exotic multi- 
pole moment in the solid state physics. 

Although it has been very rare to observe phenomena 
related to octupole moment in ri-electron compounds, 
we have sometimes found octupole-related phenomena 
in /-electron compounds. For instance, in order to un- 
derstand anomalous properties of CeBg, significant roles 
of octupole moment have been discussed actively. In par- 
ticular, apparent contradiction between neutron scatter- 
ing and NMR measurements on the so-called phase II 
of CeBg has been clearly resolved by considering effects 
of magnetic-field induced octupole moment on the hy- 
pcrfine coupling. 3 In recent years, furthermore, possibil- 
ity of ordering of octupoles themselves has been inten- 
sively discussed for CexLai-zBg 4-9 and for NpO2. 10 ~ 17 
In both materials, the crystalline electric field (CEF) 
ground states are r 8 quartets, 18-22 and in such CEF 
ground states, there are dipole, quadrupole, and octupole 
degrees of freedom. 23 In fact, several experimental facts 



have been reconciled by assuming F$ u antiferro-octupole 
ordering in these materials. 

In order to contribute to the further development of 
multipole physics, there are two different ways for the 
direction of research. One is to clarify the effect of oc- 
tupole moment in a paramagnetic phase, not in the or- 
dered state, in a phenomenological level. It is one of nat- 
ural extensions of the research for mutipoles, since the 
phenomenological theory for octupole ordering has been 
developed so far. Here it should be noted that octupole 
moments are magnetic, that is, parity of octupoles is odd 
under time reversal. Thus, in principle, octupoles would 
influence magnetic properties even in the paramagnetic 
phase. However, such effects have not been studied se- 
riously. We should make efforts to include such effects 
to analyze magnetic properties of materials with active 
octupole degrees of freedom. 

Another direction of the research is to develop a mi- 
croscopic theory for multipole-related phenomena. The 
construction of such a microscopic theory was consid- 
ered to be difficult, since it has been believed that it is 
standard to use the LS coupling scheme to express the 
/-electron state. However, we have recently developed a 
microscopic theory for /-electron systems by exploiting 
a j-j coupling scheme. 24-26 The ground and low-energy 
excited states can be well described by a microscopic 
model on the basis of the j-j coupling scheme, in contrast 
to the naive belief for the limitation of the j-j coupling 
scheme. For instance, we have successfully explained the 
microscopic origin of octupole ordering of Np02 even in 
a simple d-electron-like approximation. 25 ' 

Since we have shown that the microscopic model on the 
basis of the j-j coupling scheme is useful to understand 
multipole ordering for /-electron systems, it is natural 
to investigate further the effect of multipole fluctuations 
in a microscopic level. We believe that it is challeng- 
ing and important development of the microscopic the- 
ory for /-electron materials. In fact, it has been actually 
proposed that exotic superconductivity in PrOs4Sbi2 27 
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may be induced by quadrupole fluctuations, 28-30 
a quadrupole ordered phase appears under a magnetic 
field. 31,32 In order to discuss such a possibility from a mi- 
croscopic view point, it is necessary to include multipolc 
fluctuations in a microscopic theory for superconductiv- 
ity. In particular, we should classify complex fluctuations 
in /-electron systems into multiple components. 

In this paper, we exhibit our new development for mul- 
tipole physics including review of our previous work. In 
§ 2, first we briefly review our simple phenomenological 
theory to understand the effect of octupole moment in a 
paramagnetic phase. 33 In particular, we discuss effects of 
octupoles on the magnetic susceptibility, to clarify how 
octupoles influence macroscopic properties in the param- 
agnetic phase. In § 3, we explain the recent development 
of our microscopic theory for multipole phenomena. Af- 
ter a short explanation of the model on the basis of the 
j-j coupling scheme, we briefly review the microscopic 
theory for multipole ordering. 25 ' 26 In particular, we find 
the r 5 „ longitudinal antiferro-octupole ordering for an 
fee lattice, which has been proposed for NpGv Then, we 
explain new development of our theory to treat multipole 
fluctuations. We present our framework for the micro- 
scopic approach to evaluate multipole response functions 
with some calculated results. Finally, in § 4, we summa- 
rize this paper. Throughout this paper, we use such units 
as fi,=fcB=l- 

2. Phenomenological Approach 

In order to clarify complex effects of multipole mo- 
ment, it is useful to consider first the problem in a 
phenomenological level. In fact, there are lots of litera- 
tures treating multipole ordering along this direction and 
the phase diagram including multipole ordered state has 
been successfully explained by applying a mean-field ap- 
proximation to the Heisenberg-like model for multipole 
moments. In this section, we briefly review our recent re- 
sult based on such a phenomenological model, to clarify 
significant effect of octupole moment on magnetic prop- 
erties in the paramagnetic phase within the mean-field 
approximation. 33 

Let us here consider a cubic system for simplicity. 
There is no essential difficulty to extend the theory to 
other cases. In a cubic system, the octupole moment T 
with the same symmetry Ti u as the dipole moment J 
is an independent degree of freedom in a certain CEF 
ground state, e.g., Tg, ground state. In such a case, oc- 
tupole and dipole moments interact each other. In gen- 
eral, we should note that T is not orthogonal to J. Thus, 
we redefine T' , a linear combination of T and J, which 
is orthogonal to J, i.e., tr(J Q T^)=0, where a indicates 
a cartesian component and tr(- • • ) denotes the sum of 
the expectation values over the CEF ground states. We 
normalize T' so as to satisfy tr( J a J a ) = tr{T' a T' a ). 

Then, a mean-field Hamiltonian under a magnetic field 
H along the z-direction up to 0(H) is given by 

Hmy = - HM Z - (g Jf i B ) 2 [A d J Z (J Z ) + XoT' z (T' z ) 

+ \ do (J z (T' z )+T' z (J z ))}, 

where M z denotes magnetic moment, gj is the Lande's 



400 



(2.1) 




100 150 

r(K) 



300 



Fig. 1. Temperature dependence of the inverse of the magnetic 
susceptibility \ of U0 2 under H=V) kOe and H || [100] 36 (solid 
curve). The dotted curve represents the Curie- Weiss law with 
ju efl = 3.23^ B [C = fj,l H /(3k B )} and 9 = -223 K. The inset 
shows a schematic view of the dipole moment J and octupole 
moment T in a Ts CEF ground state. 
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Fig. 2. Temperature dependence of the inverse of the magnetic 
susceptibility \ of Np0 2 under H=10 kOe and H || [100] 39 
(solid curve). The dashed curve represents the fit obtained with 
cq. (2.5) with fi cB = 3.20^b, 6 = -146.1 K, 6 1 = -161.9 K, and 
02 = —22.5 K. The dotted curve represents the Curie- Weiss fit 
with fi c f{ = 3.20/^b an d 9 = —139 K. The inset shows a schematic 
view of the dipole moment J, octupole moment T, and T' which 
is orthogonal to J in a Fg CEF ground state. 



g-factor, and /ib is the Bohr magneton. Phenomenologi- 
cal dipolc-dipole, octupole-octupole, and dipole-octupole 
interactions are expressed by Ad, A , and Ado- We have 
multiplied the factor (g,//iB) 2 to the interaction terms 
for convenience. In general, there are other multipole in- 
teractions, 34 but the mean field from them cancel out 
within O(H) due to the cubic symmetry. Thus, eq. (2.1) 
does not change, even if we consider other multipole in- 
teractions. 

Now we diagonalize the mean-field Hamiltonian by in- 
troducing the linear combinations of J z and T z : 

M { p = -3,/Mb cos 0(cos 0J Z + sin^), (2.2) 
Mf> = -gjfiB sin (f)(sva(j>J z - cos (j)T' z ). (2.3) 
We have normalized and so that the sum of 
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them is the magnetic moment M z , i.e., M Z L> + M z z> = 
M z = —gjH&J z . By choosing the parameter <f>, we can 
diagonalize the mean-field Hamiltonian as 

Hmf = - HM^ - MM^iM^) 

(2.4) 

-hmP -\ 2 mW(mP). 

Then, the moment MP (ikfj 2) ) follows the Curie- Weiss 
law with the Curie constant C\ = C cos 2 cf) (C2 — 
Csin 2 (j)) and the Weiss temperature 8\ = C\\\ {8 2 = 
C2A2), where C=C\ + C2 is the Curie constant for a free 
ion. The following relations for the Weiss temperatures 
without the parameter <f> may be useful: 8\ + $2=C(Ad + 
A ) and 6*i#2=C 2 (AdA — A^ ). Then, the magnetic sus- 
ceptibility is given by 

_ C\ C*2 

X ~ T-8 l T -8 2 

(2.5) 

{T — Oi)(T — 62) ' 

where 8={C\8i + C^^/C^CAd is the Weiss tempera- 
ture, i.e., the magnetic susceptibility follows the Curie- 
Weiss law x—C/(T— 8) asymptotically at high temper- 
atures. 

As a typical application of the present theory, we com- 
pare experimental results of UO2 35 ' 36 and NpC>2. 37 ~ 39 In 
U0 2 , the CEF ground state is a T 5 triplet. 40 - 41 In this 
CEF ground state, the r 4 „ octupole moment is propor- 
tional to the dipole moment (see the inset of Fig. 1), and 
thus, we expect Curie- Weiss behavior. In Fig. 1, we show 
the temperature dependence of the magnetic susceptibil- 
ity of UO2 , which follows the Curie- Weiss law very well 
above the antiferromagnetic transition temperature 30 K 
at least up to 300 K. 

In Np02, the CEF ground state is a r 8 quartet. 21, 22 
In this CEF ground state, the T 4u octupole moment is 
linearly independent of the dipole moment (see the inset 
of Fig. 2), and the present theory is applicable. Indeed, 
as shown in Fig. 2, the magnetic susceptibility deviates 
from the Curie- Weiss law below 150 K. This temperature 
is much higher than the transition temperature 25 K. 
Moreover, this transition is probably a transition to an 
ordered state of T 5u octupole moments, 12-17 which do 
not influence the magnetic susceptibility at least within 
the mean-field theory. Thus, this deviation of the mag- 
netic susceptibility from the Curie- Weiss law is unlikely 
to stem from critical phenomena. Rather, the present 
theory can naturally explain this deviation, as shown in 
Fig. 2. 

3. Microscopic Approach 

In the previous section, we have briefly reviewed the 
phenomcnological approach to multipole physics, by fo- 
cusing on the effect of octupole moment on magnetic sus- 
ceptibility in the paramagnetic phase. In order to obtain 
a deeper insight into the phenomena concerning multi- 
pole degrees of freedom, it is highly requested to pro- 
ceed to a microscopic theory beyond the phenomenolog- 
ical level. For this purpose, we have recently developed 
the theory on the basis of a j-j coupling scheme. After 
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the explanation of the starting model, we show our mi- 
croscopic theory for multipole ordering as well as recent 
development for the treatment of multipole fluctuations. 

3.1 Model in a j-j coupling scheme 

In order to include the effect of Coulomb interactions, 
spin-orbit coupling, and CEF potential in /-electron sys- 
tems, the LS coupling scheme has been frequently used. 
However, it is not possible to apply standard quantum- 
field theoretical technique in the LS coupling scheme, 
since Wick's theorem does not hold. In order to over- 
come this difficulty, it has been proposed to construct 
a microscopic model for /-electron systems by exploit- 
ing the j-j coupling scheme, 2 where we include first the 
spin-orbit coupling so as to define the state labeled by 
the total angular momentum j. For / orbitals with an- 
gular momentum £=3, we immediately obtain an octet 
with j=7/2(=3+l/2) and a sextet with j=5/2(=3- 1/2), 
which are well separated by the spin-orbit interaction. 
Since the spin-orbital coupling is, at least, in the order 
of 0.1 cV for / electrons, it would be enough to take 
into account the j=5/2 sextet, when we investigate low- 
temperature properties of /-electron compounds in the 
j-j coupling scheme. 

Then, the Hamiltonian is written as 

r.a.ii.u v.a.u 

(3.1) 

+ (1/2) -^Ml,M2,A'3,^4 a r^i a r^2 aT 'A'3 a '-Ai4 7 

where a rfl is the annihilation operator of the / electron 
with the z-component n of the total angular momentum 
j=5/2 at site r, a denotes a vector connecting nearest- 
neighbor sites, t° is the hopping integral between elec- 
tron for v-st&te at site r + a and that for /n-state at r, 
B^ v is the CEF potential, and the Coulomb integral I 
in the j-j coupling scheme is expressed by using three 
Racah parameters, E , E\, and E 2 - In this study, we use 
a simple form for the hopping integrals, i.e., we consider 
the hopping integrals through the (ffcr) bonding. 

3.2 Multipole ordering 

Concerning multipole ordering, we have already re- 
ported a microscopic approach by considering correct 
/-electron symmetry. 25, 26 Thus, in this subsection, we 
briefly review the method and obtained multipole or- 
dered states. 

Let us consider again a cubic system for simplicity. 
In this case, it is convenient to express /-electron state 
in terms of the eigenstates, Tj and Tg, in a cubic CEF 
potential. Annihilation operators for the Tg states are 
given by 

/r«T = VW^ a r5/2 + vV6 a r--3/2, (3.2a) 

frai = V~5/6a r _ 5 /2 + \/T/6a r3 /2, (3.2b) 

and 

/r/3T = a rl/2, (3.3a) 

frfJi = a r -i/i, (3.3b) 
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where we introduce pseudo-spin index u{=] and {) to 
distinguish the states in each Kramers doublet, while or- 
bital index r distinguishes degenerate Kramers doublets. 
Annihilation operators for the states are given by 

(3.4a) 



/r 7 t = V l/6flr5/2 ~ V^>/^ a r-3/2 



l/6a r 



-5/2 



5 /6a 



r3/2- 



(3.4b) 



Now we consider the reduction of the Hamiltonian by 
suppressing orbital, in particular, for actinidc diox- 
ides. The reasons are as follows: First, wc note that the 
r> wave- function extends along the [111] direction and 
other equivalent directions. Since in actinide dioxides, 
oxygen anions locate in these directions, the T-j state 
should be energetically penalized. Thus, it is expected to 
be higher than the r 8 level. Second, when we accommo- 
date two, three, and four electrons in the r 8 orbitals, the 

(2) 

ground states are T 5 , T y s , and Ti, respectively. In fact, 
these states have been found as the CEF ground states 
in U0 2 , 40 ' 41 Np0 2 , 21 ' 22 and Pu0 2 , 42 ' 43 respectively. Fi- 
nally, the CEF excitation energy in Pu0 2 is as large as 
123 meV. It is expected to be the same order of magni- 
tude as the splitting between the r 8 and T7 levels, and 
thus we can ignore the T7 states at low temperatures. 

After transforming the Hamiltonian cq. (3.1) to that 
in the r7-r 8 basis, we simply discard the terms concern- 
ing orbital. Then, we obtain the r 8 -orbital Hubbard 
model as 



E 



\t'<t' frTafr+ar'o-' 



r.a.r.cr.T' 



+ ^E 



^"T"7" "J" ^TT J, 



(3.5) 



1 Srjio 



+ J ^ ' fl T pfl T ^.frT' lfrr'1 1 

where , , is the hopping integral of an electron 
with (t',0 7 ) at site r+a to the (t, <t) state at r, n rTa 
=frrafrTa, and n rT = J2<j n rr<r- The coupling con- 
stants U, U', J, and J' denote the intra-orbital, inter- 
orbital, exchange, and pair-hopping interactions, respec- 
tively. These are expressed in terms of Racah parameters, 
and we obtain the relation U=U' + J + J' . 

In order to estimate interactions among multipoles 
from the above r 8 -orbital model, we apply the second- 
order perturbation theory with respect to the hopping 
integrals of r 8 electrons, in a manner to estimate the 
superexchange interactions in (i-electron systems. Here 
we note that the hopping integrals are depend on the 
pseudo-spin a, orbital r, and the direction of the hopping 
a due to the /-electron symmetry. Thus, the hopping in- 
tegrals depend on the lattice structures, and we obtain 
different multipole interactions for different lattices. 

The obtained multipole interaction models have been 
further analyzed by using the mean-field theory. For a 
simple cubic lattice, we find a T 3g antiferro-quadrupole 
ordering at a finite temperature, and as we lower the tem- 
perature further, we find another transition to a ferro- 



magnetic state. For a bec lattice, we find a T 2u antifcrro- 
octupole transition, and a ferromagnetic transition fol- 
lows it. For an fee lattice, we analyze the multipole in- 
teraction model carefully, since the fee lattice has ge- 
ometrical frustration. Thus, we first evaluate multipole 
correlation functions numerically, and determine relevant 
interactions to the ground state. Then, we apply mean- 
field theory to a simplified model including only these rel- 
evant interactions, and find a T 5u longitudinal antiferro- 
octupole ordering. This T 5u ordering has been proposed 
for the ordered state of Np0 2 , in which Np ions forms 
an fee lattice and the CEF ground state is a r 8 quartet. 

3. 3 Multipole fluctuations 

In the previous subsection, we have reviewed our mi- 
croscopic theory for multipole ordering. It is considered 
to be natural to attempt to extend our theory to treat 
multipole fluctuations. In order to clarify multipole fluc- 
tuations from a microscopic view point, again we resort 
to a method which has been widely used in the d-electron 
research field. Namely, we apply the fluctuation exchange 
(FLEX) approximation to the model based on the j-j 
coupling scheme. 

Although there is no essential difficulty to consider 
a three-dimensional lattice, here we consider a two- 
dimensional square lattice for simplicity. By keeping both 
T 7 and r 8 orbitals, we rewrite the model Hamiltonian 
(3.1) as 



^ — 5^ VffjT'ir'/rTff/' 



r+ar'a' 



r.a.T.<T,T' 



+ ^E 



n rTl n rTl + U'/2 ^2 



Tlj"j-TL>p 



(3.6) 



where U=U'=E . Note that for simplicity, we consider 
the situation in which one of Racah parameters Eq is 
finite and the other Racah parameters are zero, since the 
interactions between them are complicated in general. 
We believe that characteristics of multipole fluctuations 
will be grasped even in this simplified model, as long as 
we treat the /-electron symmetry correctly. In eq. (3.6), 
we have introduced CEF parameters A and A', where A 
denotes the splitting energy between T 8 orbitals, while A' 
indicates the splitting between Tg and orbitals. Note 
also that A is introduced to include partly the effect of 
a tetragonal CEF potential. 

On a square lattice, the hopping integrals of / electrons 
through (ffa) bonding do not depend on the pseudo- 
spin a, and they are given by 



and 



t (a,0) 



t (0,a) 



-1/V3 






(3.7) 



(3.8) 
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where t = 9(//er) /28 and a is the lattice constant. In this 
subsection, the energy unit is t and a is taken as unity. 

In the present model, since the pseudo-spin a is a con- 
served quantity, the Green's function does not depend 
on a. Then, the Green's function is represented by 3 x 3 
matrix. The Dyson-Gorkov equation in a matrix form is 
given by 

G{k) = G (0) (fc) + G (0) (fc)S(fc)G(fc), (3.9) 

where G(k) is the Green's function and G^°\k) is the 
non-interacting Green's function. Here we introduced 
the abbreviation k=(k,ie n ), where k is momentum and 
e n =(2n + l)nT is the Matsubara frequency for fermions 
with an integer n and a temperature T. The self-energy 
E(fc) is given by 

5W*0 = § E Grlri(k-q)V Tiri;T2Ti (q), (3.10) 

where N is the number of lattice sites, q=(q,iw m ), and 
Lo m = 2mirT is the Matsubara frequency for bosons. The 
fluctuation exchange interaction V(q) is given by 

V(q) = lv°(q)+ 1 -V c (q), (3.11) 

where 

V s (q) = U s X s (q)U s - \u s X (0) {q)U a + U s , (3.12) 

V c (q) = U c X c (q)U c - l -U c x m {q)U c - U c . (3.13) 

The matrices U s and U c are given by U^ T . TT = U!f T . TT = 
U, U S TT ,. TT , = -U C TT , :TT , = U', and U^. T , T , = 2U>, where 
r 7^ t', and the other matrix elements are zero. The 
susceptibilities x s (?) an d X c (<z) are given by 

X s (q)=X {0) (q)[l-U s X {0) (q)}- 1 , (3.14) 

X c (q)=X i °\q)[l + U c x {0) (q)}- 1 , (3.15) 

where 

X ( X,r 3 rM = -^E G ^ 3 (H 9 )G V2 (fc). (3.16) 

k 

We solve the Dyson-Gorkov equation (3.9) numerically. 

When we obtain the self-consistent solution, we can 
also evaluate the multipole response functions. The mul- 
tipole response functions are given by 

X T -< (9, = £ [ dre-^^ (T T O r r -< (r)O^ (0)) , 
r Jo 

(3.17) 

where denotes the origin and j3 = l/T. The Hciscnberg 
representation for the operator O r ^ is defined as 

T r -'{r) = e m-^) T, e - m -^N)^ (3 lg) 

where \x is the chemical potential and TV = ^ r ThfT 1& 
the number operator of electrons. The multipole operator 
at site r is given by 

Or' = E frT 1 a 1 O i r^ ai - T2( j 2 frT2<y2 ■ (3.19) 

TiT2(7l(72 



The 6x6 matrix O t ~< is obtained by using the operator 
equivalent method for each multipole operator. We nor- 
malize O v i so that the sum of squares of eigenvalues is 
unity, for convenience. It is useful to introduce the 3x3 
matrices O sT ~' and O cT ~< which describe spin and charge 
parts of the matrix O v ~t , respectively: 

°rA = 4 E ^ i;T2CT2 <T CT2CTl) (3.20) 

<Ti<T 2 

Ot!T 2 = ~~tk E < ^r7cri;r 2 CT 2 ^o- 2 CTn (3-21) 

V <Ti<T 2 

where <x are the Pauli matrices. In the FLEX approxi- 
mation without the vertex corrections, the multipole re- 
sponse functions are given by 

x {q) ^ ] ^tit 2 ;t3T4 (q)o T2T1 • o T3T4 

rir 2 r 3 r4 (3^2) 
+ Y C (o)0 Cr ^ O cT -< . 

AtiT 2 ;T 3 T 4 Vyy^T 2 Ti ^T 3 T4 

In the following, we show results for a 32 x 32 lat- 
tice, at U=U'=At, T=0.02t, and A'=0. The number of / 
electrons per site is one. In the calculation, we use 1024 
Matsubara frequencies. In Figs. 3 and 4, we show the 
multipole response functions for A=0 and A=2t, respec- 
tively. Even for A'=0, the number n ri of electrons in 
the F7 orbital is almost zero for these model parameters, 
since the T-? orbital in the present model on a square 
lattice is localized and the system cannot gain kinetic 
energy by increasing the number of electrons in the 
orbital. Thus, for these parameter sets, the present model 
is virtually equivalent to the r 8 model, which has been 
already studied by Takimoto et al. 44 The antiferromag- 
netic fluctuations pointed out by Takimoto et al. can be 
classified into multipole fluctuations in /-electron sys- 
tems as shown in Figs. 3 and 4. 

For A=0, there is no significant structures in the mul- 
tipole response functions, as shown in Figs. 3 (b)-(d). 
By changing the CEF parameter A, the Fermi surface 
changes [see Figs. 3 (a) and 4 (a)], and the fluctua- 
tions around the M point [q — {tt,7t)] become signifi- 
cant. Then, the multipole response functions for the Y± u 
dipole, T^u octupole, and T§ u octupole moments enhance 
around the M point. On the other hand, the other mul- 
tipole response functions do not enhance significantly. 
The dipole and octupole components around at M 
point are naturally considered to compose antiferromag- 
netic fluctuations, but the enhancement of T 5u octupole 
fluctuation is quite peculiar. It seems to be interesting 
to seek for anomalous component of Cooper paring in- 
duced by r 5u octupole fluctuation. This is one of future 
problems. 

4. Summary 

In this paper, we have discussed three topics on mul- 
tipolcs in /-electron systems. First, we have explained 
effects of octupole moments on the magnetic susceptibil- 
ity in the paramagnetic phase of a cubic system. When 
the T^u octupole moment is independent of the dipole 
moment as in the r 8 CEF ground state, there are inter- 
actions between the octupole and dipole moments in 
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general. Then, the magnetic susceptibility deviates from 
the Curie- Weiss law even within the mean-field theory. 
By considering such effects of octupoles, we can natu- 
rally understand the difference between the temperature 
dependence of the magnetic susceptibilities of UO2 and 
N P 2 . 

Second, we have derived multipole interactions from 
the microscopic model based on the j-j coupling scheme. 
Then, we have found several multipole ordered states, 
depending on the type of lattice structure. In particular, 
for an fee lattice, we find the longitudinal triplc-q 
antiferro-octupole ordering, which has been proposed for 
NpC>2 phenomenologically. 

Finally, we have evaluated multipole response func- 
tions by applying the FLEX approximations to the 
j=5/2 model including both the r 8 and T7 orbitals. In 
/-electron systems, the complex fluctuations, which are 
composed of spin, charge, and orbital degrees of free- 
dom, can be classified into multipole fluctuations. We 
will study such multipole fluctuations by changing the 
CEF parameters, dimensionality, and Coulomb interac- 
tions, in order to clarify superconductivity relating to the 



multipole fluctuations in the near future. 
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